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SUMMARY 

A fli^t instrument that uses electric means for measuring the 
droplet-size distribution in above -free zing clouds has been devised and 
given preliminary evaluation , in fli^t. An electric charge is placed 
on the droplets and they are separated aerodynamically according to 
their mass. Because the charge placed on the droplets is a. function of 
the droplet size, the size spectrum can be determined by measurement of 
the charge deposited on cylinders of several different sizes placed to 
intercept the charged droplets. An expression' for the rate of charge 
acquisition by a water droplet in a field of coronal discharge is 
derived. The results obtained in flight with an instrument based on the 
method described indicate that continuous records of droplet -size spec- 
tr\mi variations in clouds can be obtained. The experimental instrimient 
was used to evaluate the method and was not refined to the extent 
necessary for obtaining conclusive meteorological data. 

The desirable features of an instrument based on the method 
described are (l) The instrument can be used in clouds with tempera- 
tures above freezing; (2) the size and the shape of the cylinders do 
not change during the exposure time; (3) the readings are instantaneous 
add continuous; (4) the available sensitivity permits the study of 
variations in cloud structures of less than 200 feet in extent. 


INTRODUCTION 

In the problem of all-weather flying, the physics of clouds is of 
great importance because clouds can cause aircraft icing, interference 
with radar, turbulence, and limitations of visibility. The necessity 
for a knowledge of the physical characteristics of the icing cloud is 
evident when the fundamental variables that determine the areas of 
impingement of freezing cloud droplets on aircraft parts are examined. 

The rotating-multicylinder method described in references 1 and 2 
has been used to measure liquid-water content and droplet-size distri- 
bution in icing clouds. The operation of this method is based on the 
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principle that when cylinders of various diameters are moved through a 
cloud, the amount of water intercepted per unit projected area of each 
cylinder is dependent on the diameter of the cylinder and the inertia 
of the cloud droplets. Those supercooled droplets that strike the 
cylinders freeze onto the cylinders. The collection efficiency of 
cylinders, defined as the ratio of the amoimt of water intercepted to 
the amount of water originally containe'd in the volume swept out by the 
cylinders, has been calculated (reference 1) and is a known function of 
drop diameter, cylinder diameter, airspeed, temperature, and pressure. 
Collection efficiency is an increasing function with increasing droplet 
diameter, decreasing cylinder diameter, and increasing airspeed. From 
these known relations, the liquid-water content, the volume-median drop 
diameter^, and the drop size distribution can be calculated from the 
weights of ice collected dxiring the simultaneous exposure of several 
cylinders of different diameters. 

Data obtained by the rotating -multicylinder method have been valu- 
able in computations for the design of thermal ice -prevention systems. 
The cylinders collect ice in a manner similar to that of airplane 
components; therefore, the intrinsic value of the method should not be 
underestimated. There are limitations, however, to the rotating- 
multicylinder method in the study of cloud physics. An important 
phase in the study of the formation of clouds is the investigation 
of the drop size and drop-size distribution in clouds above the 
freezing temperature. A knowledge of the conditions in clouds above 
the freezing temperature is also necessary for the study of turbulence 
and atmospheric forecasting. 

A method that permits the study of the continuous variations in 
the structure of clouds having temperatures above freezing has been 
devised and given preliminary evaluation in flight at thfe MCA Lewis 
laboratory. The method is based on the principle that an electric 
charge of known quantity can be placed on each droplet in a sample 
of. the cloud. Because the amount of electric charge placed on the 
droplets is a function of their size, their mass can be measiared 
electrically. The droplets in the cloud sample are charged by a coro- 
nal discharge taking place between two plates oriented parallel to 
an airstream. After the droplets are charged, they are separated 
according to their size. The method of separation involves the use 
of cylinders of various sizes placed in the path of the droplets. 

As with the rotating- cylinder method, collection efficiencies of the 
cylinders of different sizes are functions of drop size and cylinder 


The amount of water in all the drops of a diameter greater than the 
volume-median drop diameter is equal to the amount of water in all the 
drops of smaller diameter. Volume-median drop diameter is often 
referred to as "mean- effective drop diameter." 
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diameter. The charge deposited on each cylinder represents a measure of 
the total nximber of droplets intercepted. The rate of water intercep- 
tion by the cylinders is measured in terms of an electric ciirrent 
rather than in the form of deposited ice, as in the case of the 
rotating-cylinder method. 

An instrument employing the method described was constructed and 
attached to the under side of a light-bomber -type airplane. Records 
were obtained of the droplet-size variations in flights through various 
types of cloud in order to prove the method feasible. All the clouds 
were above the normal freezing temperature. 

A part of the method presented in this report has appeared in a 
dissertation titled "Some Aspects of Aircraft Icing Problem and an 
Instrument for Its Study" by Rinaldo J. Brun presented to the faculty 
of the Yale School of Engineering in June 1949 in candidacy for the 
degree of Doctor of Engineering. 


SYMBOLS 

The following symbols Are used herein. (All electrical units are 
in electrostatic units in the cgs system.) 

A,B,C arbitreiry constants to be evaluated according to boouidary 
conditions 

a^ VO Ivone -median droplet radius, centimeters 

Eq electrical collection efficiency of cylinders 

Ej^ mass collection ^efficiency of cylinders 

E^ radial component of electric field strength 

Ex x-component of electric field strength 

Eq undisturbed electric field strength 

Fp radial force of attraction between electron and droplet, dynes 

j current density in ionized field 

K droplet -inertia parameter 

k electron ion mobility 
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significant dimension of obstacle, centimeters ' 
distance between anode and needle points, centimeters 
mass of droplet, grams * 

total niunber of droplets in volme swept out by cylinder 
number of droplets of group size i that Impinge on cylinder 
charge on droplet 

{ 

saturation charge on droplet 

radial distance from center of droplet to electron, centimeters 

time interval during which droplet is being charged, seconds 

free-stream velocity, centimeters per second 

electrostatic potential at any point in field 

electrostatic potential applied to cathode 

x-component of droplet velocity, centimeters per second 

rectangular coordinates 

length of cathode, centimeters 

length of sweepout field, centimeters 

electron charge 

angle between radial force component of electron attraction to 
droplet and x-axis through center of droplet, radians 

terminating angle of approach of electrons to droplet, radians 

thickness of luminous layer 

viscosity of air, grams per second per centimeter 

distance from surface of droplet to electron, as fractional 
part of radius 
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p chaxge density in ionized field 

p density of air^ grams per cubic centimeter 

Si 

p density of water droplet, grams per cubic centimeter 

a specific inductive capacitance 

T ratio of electrical charge contained in particular group in 

question compared with total of charge contained in all groups 

$ dimensionless parameter defined by equation (Cl) 

\lr . force-field function 

03 radius of tube of force lines 


PRINCIPLE OF OPERATION 

A method for measuring cloud-droplet size is presented. An elec- 
trical charge is placed on the droplets and they are aerodynamically 
separated according to their mass. A system for charging cloud droplets 
is schematically shown in figure 1. The system consists of two parallel 
plates designated A and C. A coronal discharge exists between the 
plates along the distance indicated by Y. Cloud droplets enter the 
region of coronal discharge in the direction of the air velocity and 
acquire a charge while they are in the region of coronal discharge. 

After leaving the region of coronal discharge, the droplets follow paths 
governed by the air flow around obstacles placed beyond the parallel 
plates. Cylinders of various sizes are considered practical obstacles 
because the droplet impingement on cylinders is known, as described in 
references 1 and 2. Upon striking the cylinders the droplets leave 
their chaxge on the cylinder surface. The rate of charge deposit is 
measvired by a microammeter placed in the circuit as shown in figinre 1 
or by a recording galvanometer. The current readings from several 
different sizes of cylinders can be translated into droplet-size dis- 
tribution by methods subsequently described. 

Field of coronal discharge. - The coronal discharge between 
plates A and C (fig. 1) is induced by an electric potential produced 
with a high-voltage direct -ciirrent power supply. Needles, indicated 
by B in the figure, are inserted into the cathode plate C to aid in 
obtaining a uniform coronal discharge at a much lower potential drop . • 
than would be necessary with smooth parallel plates. With the existence 
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of a coronal discharge between the plates, a layer of luminous atmos- 
phere is formed arovind the points of the needles. The luminosity 
extends out into the field a short distance X (fig. 1) . The luminous 
layer is caused by the ionization processes talking place in that region. 
These processes result in the formation of positive ions and negative 
electrons by the collision of other electrons with the air molecules. 

The positive ions move to the cathode plate C, which is at a high neg- 
ative potential, and the electrons move to the anode plate A, which in 
the experiments is maintained at ground potential. 

The current measured from each different size of cylinder can be 
translated into droplet-size distribution only if an expression is 
known for the rate of charge acquisition by a water droplet in a field 
of coronal discharge. In order to derive the expression for the rate 
of charge acquisition, a knowledge of the field strength and the forces 
of attraction between the electrons and the droplets is necessary every- 
where in the field of coronal discharge. 

The mechanism of the electrical discharge throxigh gases is very 
well described in several textbooks (references 3 to 5) . The solution 
of the equation for the field strength is very difficult for the com- 
plete region between two plates in which ions of more than one kind 
exist. In this report, a solution to the field equation is obtained 
that is complete enough to be of use. The solution is obtained by:- 
(l) the use of physical conditions that best lend themselves to a solu- 
tion of the equations, (2) the measurement of as many unknowns as pos- 
sible, and (3) the use of reasonable simplifying assimiptions . 

The derivation of the expression for the electric field strength Eq 
in the region of coronal discharge is given in appendix A. The equa- 
tion derived is 



The assumption is made that in the region between the luminous layer 
and the anode plate the electrons moving toward the anode constitute 
the only current. Calculations employing this assvimption for the case 
of corona in coaxial cylinders are in substantial agreement with e:q)eri- 
mental observations (references 3 and 5 to 7) . For the instrument \ised 
in the tests reported herein, no serious error was made in the final 
results by ignoring the effect on Eq of the current in the space 
between the plates and 

^t 

Eq = ^ (2) 

was found to be valid (appendix A) . 
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Charging of droplet. - Equations are presented for the forces that 
govern the captinre of an electron by a spherical water droplet placed 
in an ionized field, such as the field between X and the anode in 
figure 1. The spherical droplet acquires a charge by its electrical 
attraction of the electrons. The electrons follow the lines of force 
that lead to the droplet, as indicated in figure 2. The droplet ceases 
to acquire more electrons from the field when the radial force of 
attraction between the electrons in the field and the droplet is zero 
or negative anywhere in the field for all angles of approach. 

The expression for the radial component of the force exerted on 
an electron of charge c situated at a neighboring point P a distance 
r from 0 is derived in appendix B. A conducting spherical droplet 
(fig. 3) of radius a with center 0 and with negative charge Q is 
considered located at a coordinate point x, where the field Eq is 
determinable from equation (l) or (2). The derivation is based on the 
following assuniptions: 

1. The field strength Eq does not vary in the neighborhood of 
the droplet because of space charge; that is, the droplet is located 

in a uniform field. The influence of the space charge on Eq has been 
shown to be small (appendix A) . 

2. Other droplets that may be in the field are so far from the one 
under consideration that they do not influence the force on the electron. 
Under average conditions, clouds have been found to consist of droplets 
of 20 microns in diameter spaced approximately 100 diameters apart, 

3. Other electrons surrounding the droplet do not affect the 
radial component of the force between the electron at P and the 
droplet. 


4. The only charges attracted by the spherical droplet are elec- 
trons, The significant number of positive ions formed are in the lumi- 
nous region X and are swept out to the cathode. The mobility of any 
ion formed by molecules in the air is so much smaller than the electron 
mobility that the charging by large ions is of no consequence. 

The radial component of force of attraction between the droplet 
and the negatively charged electron is ^ 



( 3 ) 
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The angle 0 is measured counterclockwise from the x-axis to the 
radial-force component (fig. 3) . The first bracketed term in equa- 
tion (3) is the effect of a conducting spherical droplet in a uniform 
•field caused by the field distortion in the neighborhood of the droplet; 
the second term is the effect of the charge Q on the droplet; and the 
third term (bracketed) is the effect obtained by applying Coulomb’s law 
that the force of attraction between the electron and the ungrounded 
droplet is equal to that between the electron and its image. 


As the charge Q acquired by the sphere increases in magnitude^ 
the force Fp becomes zero for successively smaller values of the polar 
angle 0. Thus^ the charge can be imagined funneled through a cone for 
which the half -vertex angle 6 decreases as the amount of charge Q 
acquired by the sphere increases (fig. 2). The expression for the rate 
of charge acquisition by the droplet is derived in appendix B. The 
charge acquired as a f-unction of time was found to be 


Q = 



Eq + jjtt 


(4) 


Inasmuch as specific inductive capacitance a = 81 for water^ the 
effect of the specific inductive capacitance is practically negligible. 
The manner in which charge acquisition varies with time is shown in 
figure 4 for three values of j and for Eq = 10.92 electrostatic 
units per centimeter. For large values of time t 

Q„=3Eoa2 (5) 


Droplet size and collection efficiency. - The droplet size can be 
obtained from equation (4) when the time interval during which the drop- 
let remains in the region of coronal discharge and the charge deposited 
by the droplet upon striking the cylinder are known. The time interval 
is obtained from a knowledge of the air -stream velocity and the physical 
dimension Y (fig. l) . The size and the number of droplets per unit 
volume of cloud are measured in terms of the current to ground from the 
cylinders of different sizes placed directly in the path of the droplets 
downstream of the coronal-discharge field. The size of the droplets 
impinging on the cylinders is a function of the size of the cylinders. 
Both the large and. small droplets impinge on the very small cylinders^ 
whereas mostly the large drops impinge on the large cylinders. A 
method for resolving the rate of ice accumulation on a group of 
different-sized cylinders into effective droplet-size distribution of a 
cloud is presented in reference 1. The method is essentially retained 
herein^ but it is modified in order to translate electrical readings of 
the current from the cylinders into effective droplet-size distribution. 
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In a cloud the water drops are often not of uniform size. Five 
different droplet-size distribution patterns have been defined for con- 
venience in classification of clouds (reference l) . The table of dis- 
tribution patterns^ reproduced herein as table I, was adapted to cover 
some of the probable range encountered in natirre. The droplet-size 
ratios given in table I are the ratios of the average radius of the 
droplets in each subdivision to the radius Uq of the volume-median 
drop size. 

In order to use the work of reference 1 to translate the electri- 
cal readings into droplet-size distribution, the following expressions 
are derived: If, in a cloud of \iniform droplet size, Nj_ is the nxuiiber 

of droplets of radius a that impinge on a cylinder in a given time and 
N is the total number in the volume swept out by the cylinder in the 
same time, the electrical collection efficiency is 

_ NiQ _ % 

- T 

In this equation^ can also be considered as the ratio of the charge 

deposited on the cylinder to the total charge in the voliame swept out by 
the cylinder. From the usual definition of collection efficiency^ 

^i 

W ^c 

in the case of uniform droplet size (distribution A, table l) . 

The electrical collection efficiency for distributions B, C, D, 
and E cannot be evaluated as readily as for A and must be computed 
by a rather lengthy method presented in appendix C. The droplet-size 
distribution is obtained by a coii 5 )arison of observed data with computed 
values of corresponding information given in table II. The method for 
performing this operation is very similar to that described in refer- 
ences 1 and 8 for ice accretion on rotating cylinders and is briefly 
outlined in appendix C. 

Droplet deflection. - A knowledge of the deflection of the charging 
droplets as they pass through the coronal field is important in order to 
determine the best position of the cylinders, as will be discussed in 
the next section. The equation of motion is written with air resistance 
neglected and with the provision that the x-component of the velocity 
obtained with this assumption will be examined to determine whether it 
is large enough to contribute appreciably to the air resistance. 
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dv^ 
-dt = 


dv^ = - 
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+ Xq 


( 6 ) 

(7) 


where Xq is the value of x when the droplet enters the ionizing 
field. 


For t = 0.0083 second, which corresponds to the condition of an 
air-stream speed of 200 miles per hour and a coronal field 30 inches 
long, a droplet 10"^ centimeter in radius deviates hy 0.98 centimeter 
from a straight line in passing through the coronal field 
(Eq = 10 e.s .u./cm and j = 10^ e.s.u./sq cm). The deviation is so 
small that very few of the droplets' are deflected against the grounded 
plates. The velocity Vx is found to be very low (in the order of 
10 ft/sec, equation (6)), so that the air resistance can be neglected 
in the equation of motion. 


A knowledge of the deflection of the charging droplets as they pass 
through the coronal field is also important in order to obtain the mass 
of the largest particles in the atmosphere, particularly when an instru- 
ment is moimted on an airplane. In equation (7), can be expressed 

in terms of the mass of a particle. In order to obtain the mass of the 
largest particle in the atmosphere, the air-stream speed is decreased 
until no ciorrent is noted from the small cylinders. The mass of ‘ the 
largest particle is calculated from equation (7) by determining the 
trajectory of particles that enter the ionizing field near the cathode 
and are deflected by the electric field in the distance Y + Yg to 
just miss the cylinders (dashed line in fig.l). 


DESCRIPTION OF INSTRUMENT 

An instrimient based on the principle of operation described in the 
preceding section was constructed. The instrument comprised four air- 
stream channels back to back as shown in figure 5. Needles were 
inserted through the cathode plate to aid in obtaining a uniform coronal 
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discharge. The needles protruded through both sides of the cathode 
plates; thus each plate acted as a cathode for two air-stream channels. 
Cylinders of four different sizes were installed directly behind the 
region of coronal discharge (fig. 6). The assembly was hung from the 
under side of the fuselage of a twin-engine light -bomber -type airplane 
(fig. 7). 


In order to obtain maximum coronal discharge in the ionizing field 
with the direct -current voltage supply available, the physical dimensions 
of an instrument must be carefully considered. The distance I from 
the needle points to the grounded anode was experimentally determined 
by constructing a laboratory model in which the current density was 
observed when Z was varied. A value of 7.5 centimeters was deter- 
mined to be the largest consistent with steady coronal discharge for 
the direct -current potential of 25,000 volts available to apply across 
the plates. The adjustment of the distance Z between the cathodes 
and the ground plates was made with the aid of numerous plugs located 
to measure the current density j at various positions (fig. 6) . The 
plugs, all of the same area, were insulated from the ground plate. The 
cathodes were adjusted with respect to the ground plates so that all 
the plugs read the same value of J . 

The length Y of the coronal dis(?harge (fig. l) was a compromise 
between the amount of charge desired on a droplet and structural com- 
plications in mounting the instrument on the airplane as well as aero- 
dynamic considerations in flight. It is desirable to charge a droplet 
nearly to saturation or to a state represented by points beyond the 
knee of the charging ciorves, such as point A on figure 4; errors in 
measuring J, and particularly in the airplane speed, which determines 
the time a droplet remains in the charging field, then reflect less on 
calculated values of Q than if the charging time were small. When 
the instrximent is used on a modern fast-flying airplane, the length 
reqviired for the desired charging time complicates the construction. 

In the instriunent, a length Y of 30 inches was chosen to allow mount- 
ing on one of the bomb -bay doors of the light bomber. 

The insulation between the cathode, which was at a high negative 
potential, and the gro\md presented some difficulties caused by moisture 
collection on the surfaces of the insulators. The insulators used were 
commercial ceramic, high-voltage transformer bushings with 3/4-inch 
holes through their centers . In order to avoid, moistiare formation on 
the insulator surfaces inverted cups were placed around the insulators, 
as shown in figure 6, to protect them from droplets that might impinge 
on them directly. A blanket of dry hot air bled from the thermal de- 
icing system of the airplane was maintained inside the cups to avoid 
condensation on the insulator sinrfaces. 
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A gasoline-engine-driven generator furnished the 220-volt, three- 
phase power required to operate the high-voltage direct -current power 
supply. Both the engine -generator and the high-voltage power supply 
were moxmted inside the bomb hay of the airplane. The high-voltage 
power supply was designed and constructed to operate at mayiTnuni voltage 
up to 25,000-foot altitude without internal arcing. In the airplane, 
the power-supply box was pressurized at ram pressture with clean hot 
air bled from the thermal de-icing system. Pressurizing the box pre- 
vented any fuel fumes that might have been present in the bomb bay from 
diffusing into the box. 

The trajectories of droplets that entered the ionized field near 
the luminous region were computed with the use of equation (7) and 
used to position the cylinders with respect to the cathode. The mag- 
nitude of the luminous layer X was observed in the darkness to be 
less than 2 millimeters for the conditions of == 25,000 volts and 
I - 7.62 centimeters. The width of the cylinders was made smaller than 
the width X by an amount such that, when the cylinders were properly 
alined with respect to the cathode and the anode, droplets which had at 
any time been in the luminous region or in the boundary layer of the 
anode failed to strike the cylinders. The insulators used to mount the 
cylinders were designed to act also as end plates for the cylinders in 
order to provide more nearly two-dimensional flow over the cylinders. 

The cylinder diameters used were l/8, l/2, 1^, and 3 inches. The 
number of cylinders used of each particular size was governed by the 
total charge intercepted by each size of cylinder. A larger number of 
the smaller-sized cy;linders was required in order to obtain a measurable 
current of 2 microamperes. The cylinders of any one size group were 
wired in parallel. The spacing between cylinders was large enough to 
avoid aerodynamic interference of the cylinders with each other. 
Cylinders of the same diameter where ganged together were spaced 
14 radii from center to center. 

I 

A multiple -re cording galvanometer was used to record the varia- 
tions in current from each set of cylinders. A single-stage amplifier, 
based on the cathode -follower principle, was used in the circuit between 
the cylinders and the galvanometer. The current could also be checked 
in flight with a vacuum-tube microammeter, which was switched into any 
one circuit when the current reading was desired. ' • 

Concern arose as to whether nitrogen or oxygen ions, which may 
have been present in the coronal-discharge field, could be swept out 
before they were carried back to the cylinders by the air flow. The 
mobility constant k for an oxygen molecule is 1.49 X 10^. With a 
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field strength of 10 electrostatic units, the ions move in the 
x-direction with a velocity of 1.49 X 10^ centimeters per second. At 
an airplane speed of 200 miles per hour, the ions move in the 
y-direction with a velocity of 8.94 X 10^ centimeters per second. An 
extension of the cathode plate a distance Yg of 7.5 centimeters 
(fig. 1) beyond the last source of ions was considered siifficient for 
the sweep -out field. 


< FLIGHT RESULTS 

The results obtained with the instrument described in the pre- 
ceding section indicate that the method of charging cloud droplets can 
be used for some phases of cloud studies. The pvirpose of the flights 
described herein was not to obtain conclusive meteorological data but 
rather to evaluate the method. 

The instrument was sensitive in detecting water droplets and other 
sizable particles present in the atmosphere even on a day apparently 
clear except for the usual haze noted in the horizon. At an altitude 
of 4000 feet over Cleveland, a bluish haze was visible between the 

and the ground. The instrument indicated a sizable current 
reading on the small cylinders at a true airspeed of 200 miles per 
hour. The airplane speed was decreased until no current was noted from 
the small cylinders. The mass of the largest particle which had a 
trajectory that just failed to reach the l/8-inch cylinder was calcu- 
lated with the use of equation (7). If the assumption that the parti- 
cles are spheres with a density of unity is allowed, the particles had 
a diameter of 4 microns. A portion of the particle-size spectrum was 
obtained from the value of .the current measured from the l/8-inch 
cylinders at a true airspeed of 200 miles per hovir and from the tra- 
jectories of various-sized particles as calculated with eq^tion (7). 
The results of this type of analysis indicated that 3 X 10^ particles 
per cubic centimeter, ranging in size from 2 to 4 microns in diameter, 
were present. No attempt was made to measure particles of smaller size. 
These results check very well with measiurements made by other investi- 
gators (reference 9) on condensation nuclei present in the atmosphere 
on slopes of mountains. 

At an altitude of 7000 feet the number of particles was about one- 
half the number at 4000 feet, but the size range was approximately the 
same. A region at 7000 feet was found about 30 miles southwest of 
Cleveland where the Instrument indicated nearly clear air.. Over Akron, 
40 miles southeast of Cleveland, the results were the same as those over 
Cleveland. Over either city readings of zero cvirrent from the cylinders 
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were not obtained until an altitude of 12,000 feet was attained. A 
spiraling descent from 12,000 feet over Cleveland indicated that the 
haze was stratified in two layers, the thinnest part being at approxi- 
mately 5000 feet. The temperature varied from freezing at 12,000 feet 
to 10*^ 0 at 2000 feet. Several other flights in apparently clear air 
checked the aforementioned results as to order of magnitude of particle 
size and content. 

On a later flight on a clear and warmer day (temperatures ranged 
from 23° 0 near the ground to 17° 0 at 3500 ft, and 13° 0 at 7000 ft) 
immediately after a rain, practically zero readings were obtained from 
the cylinders at all altitudes. The flight was over Cleveland and at 
altitudes as high as 7000 feet. 

A galvanometer record of the current from three size sets of cyl- 
inders registered during a flight through strato -cumulus clouds is 
presented in figure 8. One inch on the strip corresponds to about 
2100 feet in the atmosphere. With records of the type shown, variations 
occurring in 200 feet of cloud extent can be discerned. The high- 
frequency variations on the strip are attributed to the natiiral fre- 
quency of the galvanometer elements. In the two top traces, an increase 
in current is indicated by a downward motion whereas with the 3 -inch 
cylinder record an increase in current reading is up. 

Between the lines WW and XX the average liquid-water content, 
volume-median droplet diameter, and distribution are 0.3 gram per 
cubic meter, 11 microns, and B, respectively. The large variations 
recorded to the right of line XX axe caused by a cumulus section of 
the cloud. At the peak conditions the liquid-water content, volume- 
median droplet diameter, and distribution were 0.9 gram per cubic 
meter, 8 microns, and C to D or more, respectively. The data were 
taken at 6000-foot altitude and the average air temperature was about 
8° C. The current flow from the set of l/8-inch cylinders averaged 
less than 2 microamperes in the stratocumulus clouds and varied up to 
5 microaii5)eres in cmulus clouds. 

On several occasions stratus clouds were approached from above. 

On each of these occasions the instrument began indicating liquid- 
water droplets from 1000 to 500 feet above the apparent cloud. In one 
40-mile flight at 8000-foot altitude averaging about 300 feet above a 
thin stratus cloud, the current measured from the l/8-inch cylinders 
fluctmted by as much as 400 percent. Even though the cloud layer was 
apparently well defined below, columns of moistiire rose above the 
layer. 
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One flight was made in stratocumulus clouds in the vicinity of the 
Cleveland Airport and over the south shore of Lake Erie. Observations 
from the ground before take-off indicated precipitation at the base of 
some regions of the clouds because of the fuzzy appearance of some spots 
on the xmder sides; however, no precipitation reached the ground. Ceil- 
ing was reported at slightly over 5000 feet. At 3000 -foot altitude, the 
temperature was 3° C and the current from the different sizes of cyl- 
inders was of the same order of magnitude as was noted on previous 
flights through haze. The air was rough, indicating considerable tur- 
bulence. At 6500 feet, the temperature was 0° C and precipitation was 
encomitered in clunps or sheets that varied in length from about 400 to 
1500 feet (estimated). On passing throx;igh these clumps, the cylinders 
all read as much as four times the current throiigh other parts of the 
cloud. A portion of the flight record is shown in figure 9. 

Because the instrument was designed without anti-icing equipment, 
flights were necessarily conducted at atmospheric temperatures above 
the freezing point. The instrument was flown through slight icing 
conditions as a matter of cxariosity. The only ill effect of the ice 
noted was a large decrease in ciirrent density J, with a consequent 
loss of signal at both the plugs and the galvanometer. This decrease 
was attributed to the foimiation of ice beads on the needle points. 

The resistance between the cylinders and the anode was checked 
during flights throxigh clouds that wetted the instrument to ascertain 
that the moisture did not cavise an appreciable current leakage to the 
electrical ground. The resistance was always very high, in excess of 
700 megohms. Flights were also made tinder various conditions to deter- 
mine if atmospheric electricity might be partly the cause of the current 
measured from the cylinders. With the coronal discharge at zero no 
current readings were observed from the cylinders. This test indi- 
cates that the airplane always assumed the local potentials in the 
cloud and that the measurements of the instrument were those caused 
by the additional charge placed on the droplets by the coronal dis- 
charge . 


CONCLUDING REMARKS 

The charged-droplet cloud analyzer is based on a principle that 
depends in part on the air flow past the instrument. In this respect 
the instrument fits naturally into the condition encountered when it 
is used with an airplane. The location of the instrument on the air- 
plane may be improved perhaps to obtain a more representative sample 
of the clouds. The first model was placed on the bomb-bay doors for 
convenience in installation. 

A difficulty soon realized in the interpretation of the data is 
the lack of other independent methods for checking the results. The 



16 


MCA TW 2458 


accviracy in determining the liquid-water content and droplet-size dis- 
tribution depends on: ( 1) the error in measuring true airspeed, from 
which the total charge on each droplet, as well as the rate of cloud 
interception is obtained; (2) the error in measuring electrical values, 
partic\ilarly the galvanometer readings obtained in the form of traces 
on films (figs. 8 and 9); and (3) the errors in the assumed distribu- 
tions presented in table I as compared with the natural conditions found 
in clouds. With the instrumentation on the airplane used for the tests 
reported herein, the error in measuring true airspeed is about 2.5 miles 
per hour, which results in a negligible error in determining liquid- 
water content and volume-median droplet diameter. For the tests with 
the instrument described, the errors in measuring and reading electrical 
values varied depending on the type of cloud encountered, becaiise the 
errors involved in measuring the cxirrent density j and the applied 
voltage were insignificant as compared with those involved in 

reading the film records of the galvanometer. Between the lines WW 
and XX of figure 8, the distribution actmlly varied from A to C. At 
any one point a difference of one distribution was possible when read 
by different observers. At the peak conditions of the cumulus portion 
of the cloud (fig. 8), the distribution was C or D, depending on the 
observer. Whether the distribution was C or D’ changed the liquid-water 
content by less than 0.1 gram per cubic meter (less than 10 percent of 
the original) . The errors are much larger in resolving the data where 
the rate of change of cloud conditions is great, such as at the fringes 
of a cumulus cloud. Under such conditions it is often difficult to 
distinguish an E distribution from a B distribution. An improvement in 
the instrumentation and technique of measuring the current from the 
cylinders will improve the accuracy. ' 

The differences in- resolving the data by different observers often 
lie in that the data obtained from the films do not precisely match any 
calculated distribution curve. The difficulty in matching the curves 
is caused in part by the scattering of the data and in part by the dif- 
ference between the assumed distributions presented in table I and the 
natural conditions existing in the cloud measured. The addition of 
cylinders larger than 3 inches in diameter may lessen the uncertainty 
caused by the scattering of the data. 

An instrument based on the method of charging the cloud droplets 
has the following desirable features: 

1. The method can be applied in clouds with ten^ieratures above 
freezing and perhaps in subfreezing temperatures with proper anti- 
icing equipment. 
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2. The size and the shape of the cylinders do not change during 
the exposure time as in the rotating-multicylinder method. 

3. The uncertainty caused hy the bouncing off of droplets should 
be much less than ih the rotating-multicylinder method because the 
transfer of charge from the droplets to the cylinders is nearly 
instantaneous and what happens to the droplets after the charge is 
transferred is unimportant. 

4. The readings are instantaneoixs and continuous. 

5. The available sensitivity permits the study of the local vari- 
ations in clouds. 

6. Weighing errors are eliminated. ' 

The analysis of clouds is very difficult with any method requiring 
cylinders or other obstacles to separate the droplets according to 
their mass, because both the liquid-water content and the droplet sizes 
in the clouds are iinknowns expressed implicitly in the equations 
employed. An independent method for measuring liquid-water content in 
the clouds would permit a simplification in the reduction of the data 
to droplet size distribution. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, Febniary 16, 1951. 
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APPENDIX A 

SOLUTION OF ELECTRIC -FIELD EQUATIONS 

The problem to be solved is the determination of the field strength 
at points in the field between two parallel plates A and C (fig. l) with 
a total applied potential imposed between them. The potential drop 

is large enough to cause a coronal discharge between the plates. Many 
atteupts to solve the field equations for different ion fields are 
described in reference but usually simplifications have been neces- 
sary to such a degree as to make the results inapplicable to the present 
problem. In the present report a solution to the field equations is 
obtained that is complete enough to be of use. The solution is obtained 
by; (l) the use of the physical conditions that best lend themselves to 
a solution of the equations, (2) the measurement of as many TJ n k n owns as 
possible, and (3) the use of reasonable sinplifying assumptions. 

The variation of electrical field strength with position is 
obtained only in the region between the luminous layer and the anode 
plate. The assumption is made that in this region the electrons moving 
toward the anode constitute the only current. This assumption may be 
somewhat debatable, but it is in agreement, within experimental errors, 
with probe tests. Calculations employing this assumption are in sub- 
stantial agreement with experimental observations. These calculations 
were of the type first made by J. S. Townsend (reference 6) and J. J. 
Thomson and G. P. Thomson (reference 3) for the case of corona in 
coaxial cylinders. 

In the luminous layer, the current is the sum of that due to the 
electrons moving toward the anode and that due to the positive ions 
moving toward the cathode. Although an exact solution would be desir- 
able, it has not been found necessary for solution of the problem for 
the charged droplets. An observation in darkness of the magnitude of X 
is all that is necessary. A knowledge of the magnitude of X is nec- 
essary in order to position the cylinders properly. 

If the field in the nonluminous region is assumed independent of 
the y- and z -directions, Poisson’s equation for the electric field 
becomes 


^ = -4»p (Al) 

dx 

The assumption that Sv/Sy and Sv/^z are zero is not strictly correct 
near the needle points; however^ the cathode plate aids in obtaining a 
uniform field over most of the distance from to the anode. The 
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assumption is made that the charge density p is composed only of 
negatively charged electrons. This last assumption also iE 5 )lies that 
the current density is constant in the region considered. 


Poisson's equation will first he used to solve for the field 
intensity Eq, which is here defined as the force at a point in the 
electric field that would he exerted on a unit positive charge. From 
the usual definition of the potential, the expression is obtained: 


£ 


0 - 


dx 


(A2) 


The cxirrent ' density at any point x is 


j = 





(A3) 


where 

j flow of positive electricity per unit cross-sectional area in 
unit time 

v'^ drift velocity of electrons 

Experiments have demonstrated the mobility k to he nearly a constant, 
over a considerable range of field strength, for a particular ion in a 
given atmosphere (reference 4) . 

The use of equation (A3) in Poisson's equation (Al) yields 


dx ^2 k 


(A4) 


from which by integration and substitution of equation (A2) 

,2 


/dV\ _ „ 2 _ 

- eq - 


Srtjx 2 
“T“ + ^1 


(A5) 


and 




(A6) 
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where C and are constants of integration. The constants of 

1 ^ 

integration are evaluated by applying the boundary conditions that the 
anode is grounded and that the cathode is charged to a potential - V-t. 

A question arises as to the validity of applying the boiindary condition 
that V = - V-t when x ** Z, because this condition implies that equa- 
tion (A5) also applies in the luminous region. The justification is in 
the experimental evidence (references 1 , 10, and 11, and unpublished 
work at the Lewis laboratory). For large values'‘of dv/dx, in the 
order of 2500 volts per centimeter or more, the rate of change of poten- 
tial with respect to distance is nearly the same in the luminous region 
as in the nonluminous region. The references cited present results with 
coaxial cylindrical electrodes, but there is no reason for expecting 
that ionization principles are different for the plane parallel elec- 
trodes in the problem presented herein. The error in the determination 
of Eq from applying the aforementioned boundary conditions is neg- 
ligible with high voltage drops and low values of j, as is true in 
the instrument for which these equations are derived. 

The cxarrent density j is measured with a microammeter that is 
ponnected between the ground and a known area of the anode (figs. 1 and 
6) . Several of these elements of anode area are located where plate 
and effects are least, and are insulated from the remainder of the 
anode, which is directly connected, to ground. The mobility k is 
obtained from Langevin's equation (reference 4) 

k = 0.75 U- (A7) 


where 

c/nig ratio of charge to mass of electron 
■ mean free path of electrons 

"c^ average random velocity of electrons determined from Maxwell- 
Boltzmann velocity-distribution curve 

« 

The mean free path of electrons is obtained from 

V - 4//2 L = 5.6 L 

where L, is mean free path of air molecules (reference 4) and 
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» 0.922 


5k 'T 
m ^ 


where 


k' Boltzmann constant 


T absolute temperature 

Compton's more accurate equation (reference 4) for the electron mobility 
is not used because the refinement gained thereby is unwarranted for the 
present pvirpose. 


From the application of the first bo\mdary condition that V 
for X as 0 in equation (A6) 

C — - ^ C ^ 

^2 “ 12«j ‘"I 


as 0 


(A8) 


and from the second boundary condition that v = - V-^ for x as 2 


- V 


±]s. 
t 12jt j 




(A9) 


The constant C]^ may be explicitly evaluated by expanding the term 
raised to the 3/2 power by the binomial theorem and retaining the first 
two terms of the expansion. This operation is permissible when SnjZ/k 
is small compared with Cq. The constant is 


C 



(AlO) 


The substitution of this value for Cq 


into equation (A5) results in 



( 1 ) 


which is the required field intensity. 


The effect of space charge on Eq was very small in the experi- 
mental model described in the section titled "Description of Instrument." 

1 —6 

The highest value of current density recorded was j = ^ X lO" anrpere 

= 1 X 10 electrostatic units; = 25^000 volts = 83*3 electrostatic 
units; I = 7.62 centimeters; k » 2.02 X 10® electrostatic units. For 
these values^ eq\xation (l) becomes 
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g 2 ^ (8fl)(i X 10^) X ^ 

0 6 
2.02 X 10 



from which it is seen that the term containing 
pared with the other term on the right because 
tion (l) may be written 



X is very , small com- 
0<x<l, and equa- 

\ 


( 2 ) 


For the Instrument used in the tests, no serious error was made in the 
final results by ignoring the effect of the current density on Eq; 
however, this would not be true for instruments in which j is an 
appreciable quantity. 
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APPEHDIX B 

RATE OF CHARGE ACQUISITION BY WATER DROPLET IN 
IONIZED FIELD 

The rate of acquisition of charge is derived for a droplet located 
in an ionized field. The spherical droplet acquires charge by elec-: 
trical attraction of the electrons. The electrons follow the lines of 
force that lead to the droplet (fig. 2) . In this manner, the droplet 
ceases to acquire more electrons from the field when the radial force 
of attraction between the electrons in the field and the droplet is 
zero or negative for all angles of approach. For reasons of sin5)licity, 
the theory is developed for spherical droplets that are perfect con- 
ductors and then discussed for droplets, such as water droplets, that 
have a specific inductive capacitance of a. 

A conducting droplet (fig. 3) of radius a with center 0 and 
with a charge Q is considered located at a coordinate point x, 
where the field Eq is determinable from equation (l) or (2). The 
objective is to derive the expression for the radial component of the 
force exerted on an electron of charge € situated at a neighboring 
point P a distance r from 0. The derivation is based on the assump- 
tions listed previously in the section entitled "Principle of 
Operation. " 


Derivation of Equations 

The radial component of the force of attraction between the drop- 
let and the electron is given by equation (3). For ease in evaluating 
the significance of the equation, a transformation is made in the form 
of r = a( 1+v) . In this manner the distance from the surface of the 
droplet to the electron is given as a fractional part of the radius, 
that is, by V. 




24 


MCA TN 2458 


The charges represented by € and Q are negative, and therefore 
electrons vill no longer be forced on to the droplet when it has 
acquired sufficient charge Q to make Fp zero or negative anywhere 
in the field for all polar angles 0 of electron approach. When the 
value of V is large compared to €, the last braced term of equa- 
tion (B1), the image term, contributes very little to the value for 
Fp. If a space surroxinding the droplet is considered in which the image 
terms are negligible, it becomes evident that Fp is negative for all 
values of 0 between n/2 and 3n/2 regardless of the amount of 
charge. As the charge Q acquired by the droplet increases in magni- 
tude, the force Fp becomes zero for successively smaller values of 
the polar angle 0. Thus, the charge can be Imagined fuimeled through 
a cone for which the half -vertex angle 0 decreases as the amoxint of 
charge Q acquired by the droplet increases (fig. 2) . 

The variation of Fp experienced by electrons approaching the 
droplet along the x-axis, where 0=0, for three values of Q is 
shown in figure 10. The last term in equation (Bl), which represents 
the mutual force of attraction between the electrons and the droplet, 
contributes greatly to the magnitude of Fp only for small values of 
V, that is, when the electron is very near the surface. An electron 
approaching the negatively charged droplet is acted on by a diminishing 
force of attraction, which may become zero or repulsive for suffi- 
ciently high charge on the droplet, until the mutual force of attraction 
becomes large enough to cause the attractive force to increase again. 

The droplet is considered fully charged when it has been charging long 
enough to have raised the magnitude of Q to that required to make Fp 
zero anywhere along the zero polar axis. Additional electrons will be 
attracted to the sphere only if they have sufficient kinetic energy 
from thermal agitation or previous motion to carry them through the 
region where Fp is zero. The charge thus acquired does not, in prin- 
ciple, have any definite limit j however, the number of electrons that 
have a high velocity is small and the droplet is in the region of ioni- 
zation for only a short time. This type of charging is of minor con- 
sequence in the instrument under consideration. 

The image term may be neglected without appreciable error in the 
expression for the force of attraction provided that the results apply 
to conditions in which (av)^ is much greater in magnitude than c^. 

The value of € is 4.803 X 10“^*^ electrostatic iinit (reference 12) . 

For droplets with a radius of 10 microns (10“^ cm) the error resulting 
from neglecting the image term is small for regions where V is 
greater than 10"^. Because of the mutual attraction between the elec- 
tron and the droplet, the escape of electrons from the sxirface of the 
sphere is very improbable. When the image term is ignored, the exq>res- 
sion for the electric potential at the position of the electron becomes 
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V = 



EoT cos 6 + 


Sl 

r 


(B2) 


The rate of accumulation of charge by a droplet is best determined' 
from a consideration of the function for the lines of force surrounding 
a sphere located in an electric field. From Lamb (reference 13^ pp. 127 
and 128)^ the f\mction is obtained for the lines of force, which are 
orthogonal to the equipotential lines given by equation (B2) . Lamb 
gives the force -field functions in terms of the zonal harmonics as 


-n+1 




2 n 


■^1 n+1 * ^ ' dti 

1 -n /, 2\ 

( 1 - 1 ^ ) 


dp 


and 


^3 » - (1 - ^) 

where n is the order of the zonal harmonic and p is assigned the 
value of cos 6 for a problem involving a sphere. 


The zonal harmonics (reference 13, p. 115) that suit the problem 
vinder consideration result in the following force -field function: 


i|r = ^1^1 + i|f2 + i|f3 



0 + C(l - cos 0) 


(B3) 


where A, B, and C are arbitrary constants to be evalviated by the 
boundary conditions. The x-component of the local field strength is 
obtained from the force-field function by the expression 


1 Silf 
0) 


(B4) 


where to is defined by the relation 

r^ = + y^ + z^ = 

(figs. 2 and 3). The relations 

cos 0 ss — 

r 


(B5) 


(B6) 
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and 


a 0) 

sxn e = - 


(B6a) 


are obtained from the geometry and equation (B5) . 

The first bo\mdary condition to be applied is that 




(B7) 


when r is very large compared to the, droplet radius a. The substitu- 
tion of equations (B6) and (B6a) into equation (B3) yields 


A ... . 2Bco 3Bcd^ , CxcD 

— BIX) T „ ~ a ’Z 




(B8) 


because Sr^/^co = 2co from equation (B5) . As r approaches ca, equa- 
tion (B8) reduces to 

= Aoi (B9) 

The combination of equations (B4) ^ (B7)^ and (B9) yields 

A = -Eq 

The other boundary condition is applied very near the surface of 
the sphere where r « a, x s a, and 0*0. Differentiation in equa- 
tion (B4) gives 


E - E ^ - -1- 
\ " ^0 3 ^2 

a a 


From the definition of the field strength 


Ex=(Er)feO, =3®0+4 

r=xjaa ^ 


(BIO) 


(Bll) 


because from equation (B2) 
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The first term on the right in equation (Bll) gives the effect of the 
external field Eq, as modified by the presence of the spherical drop- 
let, and the second term the effect of the accumulated charge Q. A 
comparison of equations (BIO) and (Bll) yields 

B = C = -Q 

The final form of \lr is 

ijf = -Eq sin^ 0 - Q(1 - cos 0) (B12) 

This expression is the equation for the lines of force along which the 
electrons flow. 


At some point in the field near the surface of the droplet, where 
r _ a(l+v) and V is large enoiigh to make the image terms of equa- 
tion (Bl) negligible but small enough numerically that r « a, some 
value of 0 exists for a given value of Q for which Fp is zero. 

Here the lines of force are normal to the radius; and there are no 
lines of force, which originated at some distance from the droplet, 
terminating on the droplet for 0 greater than this particular value. 
This angle is denoted by 0q, and from equation (Bl) . 

cos 0Q (B13) 

SEqS 


The equation of the line of force passing through the point where 


r ss a and 0=6, 


IS 


^Eoa 


Of 


+ Q 1 + 


Q 


= Er 


0>" 


aV' 


+ Q 


3iEQa ^ 


k - 1) 


[_ (3EQa2)^_ 

from eqiiation (B6), (B6a), and (B12) . For very large values of r and 


2 ,2 
0 ) = 3a 


1 - 




(3EQa2) J 




Q 


3EQa , 


3Eo^a 


(3Eoa^ + Q)^ (B14) 


U- -V 1 

Physically, o) in this equation is the radius at a large distance from 
the droplet of the largest force tube for which all the enclosed lines 
of force terminate oh the droplet. All lines of force outside this 
force tube do not intercept the droplet. Electrons that were initially 
enclosed by this force tube will impinge on the droplet. 
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The rate of negative charge acquisition is determined with the aid 
of equation (B14) and by solving for 




2 jn .2 


dQ 




(3Eoa2-Q)^ 3EQ^a^ 


dt 


Q = 


3na^EQjt 
Eq + jnt 


(B15) 


For large values of time t 

Qoo= 3Eoa^ (5) 

The value for Eq is obtained from equation (l). As was previously 
stated, the current density j in the instrument built was small enough 
that Eq was virtually constant in the nonluminous region. For this 
instrument, the droplets acquired the same charge regardless of the x- 
position through which they passed in the nonluminous region. If j is 
appreciable, the charge acquired depends on the value of x. In order 
that relation (B15) apply, j must be sufficiently small that the first 
assimption of uniform field in the neighborhood of the sphere is valid. 


The effect of the specific inductive capacitance a of the water 
droplet on the charge acquired can be determined by rewriting the 
potential expression (equation (B2)) as 


V 



0-1 

0+2 ^3y 


EQr 


0 

cos 0 + — 
r 


If the req'uired chaises were made in the process of development^ the 
final result would he 


^ I. . ^ o-l^ 

^ “ V ^ 0+2 J Eq +j«t 


and 
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Inasmuch as a = 81 for water, the effect of the specific inductive 
capacitance is practically negligible. 


Experimental Verifications 

An attempt was made to check experimentally the amount of charge 
acquired by metallic spheres in falling through a region of coronal 
discharge. The spheres were individually dropped through an ionized 
atmosphere that existed between cathode and anode plates, similar to the 
setup illustrated in figure 1. The spheres were caught below the 
ionized region in a small metallic cup electrically connected to a 
quadrant electrometer. With this arrangement, only spheres of 250-micron 
radius or larger acquired sufficient charge to allow accurate measure- 
ments with the electrometer. The smallest sphere used was perhaps 
20 times larger in diameter than the estimated average droplet size in 
clouds. 

The most consistent results were obtained when the spheres remained 
in the charging field long enough to acquire , nearly the satvuration 
charge, expressed by equation (5). Under these conditions the measured 
charge averaged approximately 2 percent less than the value calculated 
with the use of equation (B15) . . The time the spheres remained in the 
charging field was varied by changing the average velocity of the sphere 
through the field. The velocity was changed by varying the starting 
height above the field. If the height from which large spheres were 
dropped was so chosen that the acquired charge was equal to the charge 
acqioired by smaller spheres dropped from lower heights, the discrepancy 
between the measvired and calculated value of Q was still about 
2 percent. 

The quadrant-electrometer method used in the experiments described 
was insufficiently accurate to verify experimentally equation (B15) when 
particles as small as cloud droplets were of interest. A method employ- 
ing the principle of the mass spectrometer to separate particles accord- 
ing to the ratio of charge to mass may, perhaps, fulfill the need. In 
this spectrometer method, the charge acquired by a sphere might be 
measvired by noting the effect of a known electrostatic field on a moving 
charged sphere of known mass. 
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APPEIKDIX C 

TRANSLATION OF CURRENT MEASUREMENTS INTO DROPLET -SIZE DISTRIBUTION 

The procedure for determining the electrical collection efficiency 
Ec for the A distribution of table I is esqplained in the section enti- 
tled "Principle of Operation." For the other distributions of table I, 
the procedure is explained with the use of a sample calculation for the 
B distribution. The electrical collection efficiency E^ for the B 
distribution is essentially a modification of the mass collection effi- 
ciency calculated with the differential analyzer (reference l) for the 
A distribution. The over -all E^ of the B distribution is found from 
the sum of the contributions of each size group in the B distribution. 

In the work of reference 1 on the impingement of cloud droplets on 
cylinders moving through the cloud, three convenient dimensionless 
parameters were used: 


18 p ^LU 
$ — ^ 
W^w 



The collection efficiency is expressible in terms of the dimensionless 
quantities K$ and K, which are listed in table III. The values of 
K$ were chosen to cover a wide range of air speeds and droplet sizes. 
The exanple chosen for presentation herein is worked for a K$ value of 
200 . 

As a basis for beginning the computations, a value of 1/K = 4.0 
is assigned to those droplets in the voliome -median group size. (This 
particular value for i/K was chosen to correspond to one of the values 
of i/K arbitrarily chosen in reference 1 for the computations for Ej^. ) 
The effect of the variation of the groim size on i/k is obtained by 
dividing the assumed value of i/K (i/k = 4.0) by (a/aQ) , because 
the radius appears to the second power in equation (Cl) . The results 
are shown in the fourth column of table IV. The collection efficiency 
Ej, for each value of i/k thus computed is obtained by interpolation 
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of table III. The values of Eq corresponding to 200 are shovn 

in the fifth column of table IV . 

The next step is to determine what contribution each size group 
makes to the total Eg for the B distribution. A term t must be 
defined as the ratio of the electrical charge contained in the partic- 
ular group in question compared with the total of the charge contained 
in all the groups. For example. 



^2 ~ 


Vs 




N^a 




^ ”7^7^ 

7 “ 2 2 2 

^2^2 

where N is the number of drops in each size group and the subscripts 
refer to the group number in table IV.. The values for N are deter- 
mined from seven relations involving the percentage mass in each group 
size: 


0.05 




+ . . . Nya^' 


0.30 




Niai^ 


W7a 


yd.? 


0.05 = 


Nna 


1=^1 


Nyay ' 


Nyay 


The values of t are presented in the sixth column of table IV. 
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The E for the B distribution when K$ — 200 is the sum of the 
individual contributions of each size group: 


E = E ,T + 
c c,l 1 


+ E 

c,7' 7 


The value of Eg thus obtained is presented as the last figure in the 
seventh column of table IV and is listed in table II. The process is 
repeated for all the other values of i/k and K<1> given in table II. 


The data necessary from flight are the true airspeed U, the air 
viscosity ti, the air density and the CTirrent per unit projected 

area from each different size of cylinders. On log-log paper the cur- 
per unit of projected cylinder area is plotted as the abscissa 
and the cylinder diameter as the ordinate (fig. 11) . A standard value 
of K$ that corresponds most closely to the flight conditions is 
chosen from table II. As only an approximate value of K$ is needed 
to give a fairly accurate result, it is permissible to estimate the 
volume-median drop size so that a. first -approximation value of K$ 
can be calculated to aid in the choice of the standard value of K$. 

The five columns corresponding to five distributions of droplet sizes 
for the standard value of K$ chosen are plotted on log-log paper, Eg 
as the abscissa and i/k as the ordinate (fig. 12) . Log-log paper 
should be used because multiplication of the value on one axis corre- 
sponds merely to a displacement along that axis. Because Eg and cur- 
rent per unit of projected cylinder area are proportional and because 
1/K is proportional to the cylinder diameter 2L, the plotted points 
of figure 11 closely coincide with one of the curves of figure 12 when 
the two sheets of paper are superimposed and shifted horizontally and 
vertically (not rotated), provided K$ was estimated fairly closely. 
The best fit of the experimental curve of figiire 11 with one of the 
con 5 )uted curves of figure 12 is shown in figure 12. 


If the estimated value for K<1> is incorrect, the curve of figure 11 
cannot be made to coincide with any of the curves of figure^ 12. If the 
ciirves do not coincide, figure 11 is matched with plots of l^K against 
Eg for other values of K$. A method for interpolating K$ is pre- 
sented in references 1 and 8. 
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table I - FIVE ASSUMED DISTRIBUTIONS OF DROPLET SIZE 

[size is ejq)ressed as ratio of average drop radius 
in each group to volume-median drop radius Uq, 
Table is taken from reference ij 



Total liquid 
water in each 
group size 
(percent) 

Distributions 

A 

B 

C 

D 

E 

5 

1.00 

0.56 

0.42 

0.31 

0.23 

10 

1.00 

.72 

.61 

.52 

.44 

20 

1.00 

.84 

.■77 

.71 

.65 

30 

1.00 

1.00 

1.00 

1.00 

1.00 

20 

1.00 

1.17 

1.26 

1.37 

1.48 

10 

1.00 

1.32 

1.51 

1.74 

2.00 

5 

1.00 

1.49 

1.81 

2.22 

2.7i 


Example of interpretation: 30 percent of the liquid-water con- 

tent of any cloud is contained in droplets that have a radius Uq. 

In the case of the B distribution, 20 percent of the liquid-water 
content is contained in droplets that have a radius smaller than the 
volume-median radius Sq by a ratio a/a^ = 0.84 and another 20 per- 
cent in droplets that have a radius larger than a^ by a ratio • 
a/a^ = 1.17. A similar interpretation applies to the remaining values. 
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TABLE II - VALUES OF ELECTRICAL COLLECTION EFFICIENCY E^, 


FOR DROPLET-SIZE DISTRIBUTIONS A, B, C, D, AND E 
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TABLE III - VALUES OF Ej^ AND E^, FOR DISTRIBUTION A AS FUNCTION 


OF i/k for various values of K® 


[Table from reference ij] 


l/K 



% 

= Ec 




0 

20 

200 

1000 

3000 

10,000 

4.0 

0.040 

0.030 

0.021 

0.015 

0.011 

0.009 

2.0 

.187 

.137 

.101' 

.076 

,060 

.044 

1.0 

.385 

.326 

.271 

.226 

.186 

.149 

.5 

.560 

.522 

.478 

.420 

.373 

.314 

.2 

.761 

.729 

.694 

.656 

.624 

.582 

.1 

.864 

.843 

.822 

.796 

.767 

.736 

.05 

.927 

.916 

.904 

.886 

.870 

.848 

.02 

.970 

.963 

.961 

.951 

.943 

.932 

.01 

.984 

.983 

.983 

.985 

.971 

.965 


TABLE IV - SAMPLE CALCULATIONS FOR E^ FOR DISTRIBUTION B, 

AND K3> » 200 


1/K for VO iTime -median drop size assigned value 4,0. 


Group 

num- 

ber 

a 

^0 

Percentage 
in eacli 
size group 

1 

K 

Ec 

T 

TEc 

1 . 

0.56 

5 

12.73 

0 

0.085 

0 

2 

.72 

10 

7.72 

.003 

.133 

.0004 

3 

.84 

20 

5.68 

.008 

.227 

.0018 

4 

1.00 

30 . 

4.00 

.023 

.287 

.0066 

5 

1.17 

20 

2.92 

.050 

.163 

.0082 

6 

1.32 

10 

2.29 

.081 

.072 

.0058 

7 

1.49 

• 5 

1.82 

.118 

.032 

.0038 







Z = 0.0266 
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field. Undisturbed electric field strength Eq = 10.92 electrostatic units 
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Cathode needles 



Figure 5. - Bottom view of charged-droplet cloud analyzer. Arrows indicate direction of air stream through the channels 
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Figure 6. - Schematlo diagram of charged-droplet cloud analyzer, 
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Figure 7. - Rear view of charged- droplet cloud analyzer mounted on bomb-bay doors of airplane 
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Figure 8. - Flight record through stratocumulua clouds. Trace of fourth cylinder too faint to reproduce in this figure 
For the 1/8- inch and l/2-inoh cylindera an increase in current is measvired downward whereas for the 3- inch cylinder 
an increase in current is measured upward. 
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Figure 9. - Flight record through stratocumulus cloud congested with clumps of precipitation. Trace of fourth cylinder 
too faint to reproduce in this figure. For the l/S-inch and l/2-inch cylinders an increase in current is measured 
downward whereas for the 3- inch cylinder an increase in current is measured upward. 
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Figure 10. - Force e^qjerlenced by electrons approaching charged 
spherical droplet along x-axis for several values of Q. 


Cylinder diameter > 2L, cm 
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Figure 11. - Electric current measvired from cylinders of 
different size In cumulus cloud. Airplane speed, 200 miles 
per hour; altitude, 6000 feet. 





Reciprocal of droplet -inertia parameter, l/K 
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Electrical collection efficiency, Eq 


Figure 12. - Droplet -inertia parameter as function of electrical 
collection efficiency. K$, 200; 1 /k for voliune-median droplet, 

4.0. 
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